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( 摘要 :介绍 了 通过 Matlab/Simulink 建立 干涉 式微 波 辐 射 计 系统 仿真 模型 的 方法 。 针 对 干涉 式微 
, 波 辐射 计 的 基本 构成 单元 一 一 二 元 干涉 仪 ,搭建 了 仿真 模型 ,通过 对 辐射 计 随 机 噪声 测量 过 程 的 模 
拟 , 实 现 了 单 通道 全 功率 测量 仿真 及 双 通 道 干涉 测量 仿真 ,在 时 域 范围 内 实现 了 对 系统 测量 精度 、 


@ 灵敏 度 指 标的 定量 化 分 析 。 
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.二 对 全 球 海洋 盐 度 的 探测 工作 具有 重大 的 意义 ， 
在 助 于 人 们 进一步 深入 了 解 气候 变化 及 水 循环 等 领 
域 欧 问题 。 但 目前 关于 海洋 盐 度 的 探测 方式 还 有 很 
多 痪 限 性 , 现 有 的 研究 成 果 难 以 为 水 汽 循环 的 研究 
提 刁 有 力 的 数据 支撑 。 为 了 解决 这 一 问题 ,在 最 佳 
探测 海洋 盐 度 的 工 波 段 ,很 多 国家 都 开展 了 相关 的 
课题 研究 ,主要 包括 欧 空 局 的 探测 土壤 湿度 与 海水 
+h JF p E Æ (SMOS, Soil Moisture and Ocean Salin- 
ity)5 ,美国 测量 海水 盐 度 的 宝 瓶 座 计 划 (Aquari- 
us) 中 。 针 对 这 种 背景 ,中 国 科学 院 空间 科学 与 应 用 
研究 中 心 提出 了 主 被 动 联合 探测 式 盐 度 计 项 目 计 
A] ,在 系统 总 体 方案 上 综合 考虑 SMOS 和 Aquarius 
的 特点 ,利用 高 灵敏 度 、 高 稳定 度 的 工 波 段 一 维 干 
涉 式 综合 孔径 辐射 计 ,实现 对 全 球 海洋 盐 度 分 布 的 
高 精度 观测 ,同时 利用 共用 反射 面 天 线 的 波段 数 
字 波 束 形成 散射 计 ,实现 海面 粗糙 度 信息 的 获取 ,有 
助 于 改进 海洋 盐 度 测量 的 精度 。 
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本 文 的 研究 正 是 以 主 被 动 联合 探测 式 盐 度 计 项 
目 为 依托 ,针对 其 中 的 工 波 段 高 稳定 度 干 涉 式微 波 
辐射 计 系 统 的 基本 构成 单元 一 一 二 元 干涉 仪 , 进 行 
了 Matlab/Simulink 仿真 工作 ,搭建 出 了 仿真 模型 ， 
利用 仿真 系统 与 实际 硬件 的 一 一 对 应 关系 ,在 仿真 
FF 台 上 实现 对 系统 性 能 及 相关 指标 的 分 析 。 


2 于 涉 式微 波 辐 射 计 测量 原理 


微波 辐射 计 相 关 原 理 

根据 辐射 量子 理论 可 以 得 知 外 ,处 于 热力 学 平 
衡 的 物体 所 发 射出 的 功率 是 其 物理 温度 的 线性 函 
数 。 微 波 遥 感 的 工作 机 理 是 利用 高 灵敏 度 的 接收 机 
来 获得 目标 物质 的 辐射 噪声 功率 ,然后 通过 合理 的 
时 间 积 分 累加 ,最 终 将 辐射 功率 与 接收 机 自身 的 噪 
声 加 以 区 分 。 就 辐射 计 而 言 ,辐射 噪声 功率 的 发 
射 源 就 是 目标 物质 本 身 ,辐射 计 被 动 地 接收 天 线 收 
集 起 来 的 辐射 能 量 功率 ,可 以 通过 功率 最 终 定义 出 
物体 的 辐射 测量 亮 温 度 。 
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辐射 测量 的 灵敏 度 定义 为 辐射 计 输 出 端 可 以 
检测 到 的 天 线 辐射 亮 温 的 最 小 变化 .假设 整个 接 
收 机 系统 的 噪声 亮 温 值 为 T,, 接收 机 的 带宽 为 
也 ,系统 中 积分 器 的 积分 时 间 为 *, 则 辐射 测量 灵敏 
度 表示 为 : 


Ts 
AT = JE (1) 

从 测量 角度 来 说 ,AT 是 与 ,测量 值 有 关 的 标 
准 差 的 大 小 。 
2.2 二 元 干涉 仪 相 关 原 理 

传统 的 微波 辐射 计 在 应 用 上 有 一 些 限制 ,一 方 
面 ,高 空间 分 辨 率 的 要 求 受 限于 天 线 物理 口径 的 大 
小 , 男 一 方面 ,天 线 波束 在 每 个 分 辩 单 元 越 来 越 短 的 
FE 留 时 间 也 会 降低 辐射 亮 温 值 的 分 辩 率 '""。 基 于 这 
些 问题 ,提出 了 干涉 式微 波 成 像 技 术 。 来 源 于 射电 
天 文 领域 的 干涉 式 成 像 技 术 '" ,基本 思想 就 是 在 空 
间 频 率 域 实 现 空 间 域 的 测量 效果 号 ,将 天 线 接收 到 
的 信号 进行 复 相 关 运 算 , 输 出 目标 亮 温 分 布 的 空间 
频率 采样 点 , 称 为 可 见 度 函 数 (VE,Visibility Func- 
tion) 。 由 于 微波 辐射 计 的 最 终 目的 是 得 到 物质 的 
辐射 亮 温度 ,通过 对 可 见 度 冰 数 进行 反 傅 立 叶 变换 ， 
可 以 获得 辐射 亮 温 分 布 图 像 。 本 文 所 涉及 到 的 仿真 
模型 是 针对 上 述 过 程 中 的 第 一 步 , 即 通 过 Matlab/ 
Simulink 搭建 一 个 二 元 干涉 仪 仿真 系统 ,实现 可 见 
度 函 数 测量 过 程 的 仿真 。 

干涉 式微 波 辐 射 计 系 统 比 较 复 杂 , 目前 针对 此 
系统 的 仿真 都 是 基于 傅立叶 关系 的 空间 频率 域 和 空 
间 域 的 系统 仿真 ”"”。 然 而 ,辐射 计 测 量 的 噪声 具有 
随机 性 和 不 相干 特性 ,之 前 的 仿真 在 系统 的 性 能 分 
析 上 受到 一 定 程 度 的 限制 ,比如 对 于 量化 过 程 无 法 
仿真 ,也 不 能 进行 辐射 计 测 量 精 度 .灵敏 度 等 指标 的 
4 Wr 。 基 于 这 个 问题 ,考虑 在 时 域 上 仿真 干涉 式 
微波 辐射 计 系 统 , 通 过 在 Matlab/Simulink 平台 上 
对 二 元 干涉 仪 进行 时 域 仿 真 , 实 现 系统 功能 。 二 元 
干涉 仪 的 原理 框图 如 图 1 所 示 。 


| 输入 信号 1 | 


sr 


Li 


数字 相关 输出 


| 输入 信号 2 | 
图 1 二 元 干涉 仪 原理 框图 


Fig.1 Functional block diagram of binary interferometer 
基于 二 元 干涉 仪 的 基本 原理 ,结合 主 被 动 联合 


探测 盐 度 计 项 目 中 的 工 波段 高 稳定 度 辐 射 计 接收 
机 系统 ,利用 Matlab/Simulink 搭建 出 干涉 式微 波 
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辐射 计 仿 真 系统 模型 ,实现 可 见 度 函 数 测量 过 程 
的 仿真 。 


3 仿真 模型 


Matlab/Simulink 是 由 Mathworks 公司 推出 的 
一 种 科学 计算 仿真 软件 ,作为 一 种 可 视 化 仿真 工具 ， 
目前 已 被 广泛 应 用 于 控制 理论 和 数字 信和 号 处 理 的 复 
杂 仿 真 和 设计 5 ,前 人 经 过 探讨 实验 ,已 经 将 其 应 
用 于 仿真 模型 的 搭建 工作 中 ,并 得 到 相应 结果 。 
下 面 重 点 围绕 本 文 的 研究 方向 ,详细 说 明 仿 真 模型 
的 搭建 过 程 。 
3.1 LEER E EE HR VERE TE A ZB TT 

主 被 动 联合 探测 式 盐 度 计 中 的 辐射 计 系 统 采用 
一 维 综合 孔径 技术 , 设 定 为 8 单元 ,各 单元 接收 机 均 
为 双 极 化 "下 。 高 稳定 度 辐 射 计 接收 机 组 件 由 以 下 
几 个 部 分 构成 : 双 通 道 接收 机 前 置 .前 端 .中 频 、 公 共 
噪声 源 及 功 分 单元 .公共 本 振 及 功 分 单元 `. 数 字 相 关 
器 模块 ,系统 整体 框图 如 图 2 所 示 。 


2 工 波 段 高 稳定 度 辐射 计 接收 系统 框图 
Fig.2 Block diagram of the L band radiometer 


system with high stability 


3.2 系统 仿真 模型 概述 

接 下 来 重点 说 明 利用 Matlab/Simulink 所 搭建 
的 仿真 模型 。 仿 真 的 模型 系统 由 6 部 分 模块 构成 ， 
分 别 是 前 置 单元 、 前 端 单元 、 中 频 音 元、 公共 噪声 源 
单元 ` 公 共 本 振 单 元 以 及 相关 器 单元 ,系统 的 总 仿真 
框图 如 图 3 所 示 。 

天 线 信号 源 用 Simulink 的 白 噪 声 信 号 模块 来 
模拟 ,参数 的 设置 通过 噪声 功率 值 来 表征 相应 的 
亮 温 。 输 入 的 天 线 信号 经 过 前 置 、 前 端 、 中 频 、 相 
关 等 模块 单元 的 处 理 ,完成 二 元 干涉 仪 的 基本 功 
能 ,对 两 个 单元 天 线 信号 的 输出 进行 复 相 关 运 算 ， 
最 终 输 出 信号 的 空间 频率 采样 点 , 即 可 见 度 函 数 
VF 的 采样 值 。 接 下 来 对 系统 中 的 4 个 重点 单元 
模块 进行 描述 。 

3.2.1 前 置 单元 
前 置 单元 主要 包括 一 分 二 功 分 器 、20 dB 耦合 
器 ,仿真 框图 如 图 4 所 示 。 
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前 置 单元 1 


roa 
AWGN 


常量 2 RER 


迪克 开关 1 


图 4 前 置 单 元 仿真 框图 


Fig.4 Simulation diagram of front unit 


在 前 置 单元 中 , 主要 实现 天 线 信 号 与 公共 噪声 
源 的 混合 ,并 完成 噪声 源 注入 的 时 序 控 制 。 仿 真 在 
控制 公共 噪声 源 的 周期 性 注入 时 ,涉及 到 脉冲 时 序 
的 设 定 ,搭建 了 一 个 迪克 开关 模块 来 控制 时 序 ,具体 
来 说 就 是 通过 离散 脉冲 序列 来 实现 时 序 循环 ,迪克 
开关 的 仿真 框图 如 图 5 所 示 。 


输入 2 


切换 开关 
图 5 迪克 开关 仿真 框图 


Fig. 5 Simulation diagram of Dicke switch 


前 端 单元 1 


z 
E 


中 频 单元 1 


中 频 单元 2 
图 3 两 单元 辐射 计 系统 仿真 框图 


Fig.3 Simulation diagram of the radiometer system 


3.2.2 前 端 单 元 

前 端 单元 中 ,重点 完成 不 同 噪声 源 的 切换 性 注 
人 ,通过 迪克 开关 控制 时 序 ,将 天 线 信和 号、 公共 噪声 
源 .匹配 负载 以 及 内 部 噪声 源 接 人 接收 机 ,同时 考虑 
到 接收 机 本 身 始终 存在 噪声 ,所 以 添加 一 个 噪声 源 
模块 来 模拟 。 仿 真 中 涉及 两 个 通道 的 接收 机 ,搭建 
结构 相同 ,其 中 一 个 通道 的 仿真 框图 如 图 6 所 示 。 
3.2.3 ”中频 单元 

中 频 单 元 包括 带 通 滤波 器 、 放 大 需 、 混 频 器 以 及 
平方 律 检 波 器 等 ,主要 功能 是 实现 前 端 输出 噪声 信 
号 的 正 交 混 频 ,得 到 I/Q 两 路 信号 ,其 中 一 个 通道 
所 对 应 的 仿真 框图 如 图 7 所 示 。 

前 端 单元 输出 的 信号 进入 中 频 单 元 ,首先 通过 
带 通 滤波 器 的 处 理 实现 带 外 抑制 ,经 过 放大 后 ,与 本 
地 振荡 信号 进行 混 频 运 算 , 得 到 L/Q 两 路 信号 ,再 
经 过 中 频 滤波 、 放 大 ,输出 到 后 续 相 关 融 模块 。 男 一 
方面 ,对 主 路 输出 的 1/Q 两 路 信号 进行 15 dB 耦合 ， 
再 经 过 平方 律 检 波 以 及 低 通 滤波 器 的 处 理 , 输 出 全 
功率 信号 。 

3.2.4 数字 相关 器 单元 

数字 相关 咒 单 元 主要 包括 采样 .量化 以 及 复 相 
关 值 计算 模块 ,仿真 框图 如 图 8 所 示 。 

图 8 中 的 VF 模块 用 于 实现 数字 复 相关 值 的 计 
5t ,仿真 框图 如 图 9 所 示 。 

在 计算 复 相 关 值 的 过 程 中 ,涉及 到 信号 相 乘 、 累 
加 积分 等 数据 处 理 , 输 出 信号 的 数字 互相 关 值 ,对 数 
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常量 4 TOO 
Vries 


6 ”前端 单元 仿真 框图 


Fig.6 Simulation diagram of front-end unit 
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7 中频 单元 仿真 框图 


Fig.7 Simulation diagram of IF unit 
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图 8 数字 相关 器 单元 仿真 框图 


Fig.8 Simulation diagram of digital correlator 
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图 9 相关 仿真 框图 


Fig.9 Simulation diagram of correlation 


据 进行 后 续 校 正 处 理 ,获得 模拟 相关 系数 ,之 后 再 通 
过 可 见 度 函数 归 一 化 ,最 终 获得 空间 频率 采样 点 。 


4 仿真 参数 设 定 


4.1 系统 参数 压缩 


在 真实 的 辐射 计 硬 件 系统 中 ,接收 机 工作 在 工 


波段 ,而 要 在 计算 机 上 运行 MHz 的 时 钟 频 率 几 乎 
是 不 可 能 的 ,所 以 在 仿真 过 程 中 ,要 将 硬件 系统 
中 的 信号 频率 、 系 统 带宽 等 参数 进行 一 定 比例 的 压 
缩 。 在 确保 系统 的 采样 及 性 能 指标 满足 要 求 的 前 提 
下 ,优化 后 的 系统 参数 情况 如 表 1 所 示 。 

表 1 接收 机 参数 压缩 对 比 


Table 1 Compression parameters of the receiver 


接收 机 参数 实际 硬件 系统 仿真 模型 系统 
中 心 频率 1.413 GHz 10 kHz 

3 dB 带宽 27 MHz 2 kHz 

本 振 频 率 1.35 GHz 8 kHz 

中 频带 宽 50~77 MHz 1~3 kHz 


压缩 模型 的 原始 频段 范围 会 不 可 避免 地 降低 系 
统 的 灵敏 度 指 标 , 针 对 系统 参数 压缩 后 相 较 于 原始 
系统 的 指标 变化 ,进行 了 等 效 性 对 比分 析 。 设 定 系 
统 亮 温 值 为 280 K, 对 应 的 系统 主要 参数 指标 如 表 2 
所 示 。 


表 2 仿真 系统 与 原 系 统 指标 对 比 表 
Table 2 Comparison of the system parameters 
系统 "n 单 点 积分 观测 灵敏 度 
参数 mM 带宽 
亮 温 /K 时 间 /s 周期 /s /K 
原 系统 27 MHz 0.2 2 0. 100 
仿真 系统 2 kHz 0.5 100 0. 626 


从 2 表 中 可 以 看 出 ,仿真 中 相应 参数 的 压缩 确 
实 对 系统 的 指标 之 来 了 影响 。 其 中 辐射 计 的 灵敏 度 
指标 取决 于 系统 的 噪声 温度 、 信 号 带宽 以 及 系统 的 
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积分 时 间 。 当 系统 的 噪声 温度 一 定时 ,辐射 计 的 灵 
敏 度 与 积分 时 间 和 系统 带宽 成 反比 ,所 以 为 了 提高 
灵敏 度 指标 ,必须 增加 信号 的 带宽 值 或 提高 积分 时 
间 。 但 是 由 于 仿真 中 计算 机 性 能 的 限制 ,对 原 系统 
的 带宽 进行 了 压缩 ,不 可 避免 地 恶化 了 系统 的 灵敏 
度 指标 。 为 了 弥补 压缩 带宽 带 来 的 指标 影响 ,在 仿 
真 中 考虑 拉 长 系统 的 积分 时 间 来 优化 灵敏 度 参数 ， 
通过 在 一 个 较 长 的 观测 周期 内 对 单 点 积分 时 间 进 行 
统计 平均 ,达到 拉 长 积分 时 间 的 目的 ,有 效 地 改善 了 
系统 的 灵敏 度 指标 。 另 外 ,为 了 使 仿真 模型 能 够 产 
生 可 靠 的 噪声 亮 温 估计 ,必须 要 提高 校准 和 观察 时 
间 来 保障 长 期 稳定 的 目标 场景 温度 ,或 者 进行 多 次 
重复 试验 来 获得 平均 值 , 以 保证 准确 度 。 

4.2 提高 仿真 效率 

Simulink 的 仿真 大 多 为 非 实 时 仿真 ,对 于 一 些 
仿真 实时 性 要 求 较 高 的 场合 ,往往 存在 着 仿真 运行 
速度 较 慢 的 缺陷 , 达 不 到 预期 的 理想 效果 "5 。 为 了 
解决 这 一 矛盾 ,提出 了 利用 RTW (Real-Time 
Workshop) 实 时 代码 生成 工具 的 方法 ,将 Simulink 
模型 向 其 他 语言 模型 转换 ,以 满足 不 同 实时 仿真 速 
度 的 要 求 。 由 于 C/C 十 十 语言 的 执行 效率 较 高 , 具 
有 稳定 和 运行 速度 较 快 的 特点 ,能 够 满足 实时 仿真 
要 求 ,因此 选择 将 模型 转换 成 C/C 十 十 语言 模型 , 提 
高 仿真 运行 效率 。 

经 验证 ,通过 实时 模型 代码 得 到 的 结果 与 直接 
运行 Simulink 模型 的 结果 几乎 是 相同 的 ,表明 了 用 
RTW 生成 代码 的 正确 性 25 。 在 代码 环境 下 运行 仿 
真 模型 ,将 原始 的 仿真 速度 提升 了 几 个 数量 级 ,达到 
了 提高 仿真 运行 效率 的 目的 。 


5 仿真 结果 


5.1 单 通道 全 功率 测量 仿真 

对 于 单 通道 的 全 功率 测量 ,重点 分 析 中 频 单 元 
的 耦合 输出 信号 。 在 中 频 单 元 ,信号 混 频 后 分 成 
IQP E ,通过 平方 律 检 波 器 ,使 得 系统 的 输出 信号 
功率 值 与 信号 亮 温 值 成 线性 关系 ,再 经 过 低 通 滤 波 
器 进行 滤波 , 滤 除 信号 中 的 高 频 起 伏 分 量 , 即 输出 全 
功率 信号 。 

仿真 系统 的 参数 设置 为 :整个 系统 的 采样 频率 
设置 为 120 kHz. 仿真 时 间 210 s. 系统 中 频带 宽 
2 kHz, 中 频 耦 合 输出 端的 低 通 滤波 器 积分 时 间 
c X 0.5 s。 利 用 迪克 开关 控制 噪声 信号 的 周期 性 注 
入 时 ,采用 的 脉冲 时 序 采 样 周期 为 10 s, 每 个 时 序 周 
期 内 设 定 10 个 脉冲 值 ,通过 不 同 噪声 源 对 应 脉冲 时 
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序 值 的 设计 ,实现 系统 噪声 注 人 状态 的 切换 。 另 外 ， 
结合 实际 的 辐射 计 硬 件 系 统 , 将 系统 中 涉及 到 的 几 
种 噪声 源 信号 亮 温 值 分 别 设置 为 :天 线 信和 号 亮 温 值 
T,—80 开 , 公 共 噪 声 源 (外 品 ) 亮 温 值 T.—70 开 , 匹 
配 负载 亮 温 值 T,—300 天 ,内 部 噪声 源 亮 温 值 T, = 
100 K ,接收 机 噪声 亮 温 值 T..—200 K。 需 要 注意 
的 是 ,以 上 噪声 源 亮 温 值 均 表示 等 效 到 同一 输入 端 
以 后 的 数值 大 小 。 
5.1.1 开关 时 序 控 制 

由 于 仿真 系统 中 涉及 到 多 种 不 同 的 噪声 源 , 包 
括 天 线 信 和 号、 公共 噪声 源 ( 外 噪 ) .内 部 噪声 源 、. 匹 丁 
负载 噪声 源 以 及 接收 机 噪声 ,要 对 这 些 噪声 源 的 洽 
人 状态 进行 控制 ,就 涉及 到 多 个 迪克 开关 的 时 序 设 
置 问 题 ,不 同 噪声 源 对 应 的 控制 开关 时 序 不 同 , 下 面 
详细 说 明 脉 冲 时 序 的 设置 情况 。 首 先 , 考 虑 到 实际 
的 接收 机 通道 中 始终 存在 系统 噪声 ,所 以 接收 机 噪 
声 源 信号 始终 添加 在 系统 中 。 其 次 ,天 线 信 号 注入 
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到 系统 中 ,在 经 过 匹配 负载 噪声 源 处 的 迪克 开关 时 ， 
通过 开关 的 切换 选择 功能 ,将 系统 整体 切换 于 这 两 
种 信号 源 之 间 , 所 以 将 匹配 负载 噪声 源 处 的 控制 时 
序 设 置 为 [0000000111], 其 中 ,0? 表 示 信 和 号 不 
输入 系统 ,“1” 表 示 信 号 输入 系统 。 可 以 看 出 ,在 第 
8、9、10 这 3 个 时 序 的 情况 下 将 匹配 负载 噪声 源 注 
和 系统 ,其 他 时 序 情况 下 将 天 线 信 号 注 和 系统。 最 
后 ,在 天 线 信和 号 噪声 源 或 匹配 负载 噪声 源 这 两 种 状 
态 切换 的 基础 之 上 ,搭配 其 他 噪声 信号 的 输入 ,实现 
不 同 的 噪声 状态 ,用 以 完成 系统 定 标的 功能 。 对 于 
公共 噪声 源 ( 外 噪 ) 的 控制 开关 ,设置 的 脉冲 时 序 为 
L0000010000], 即 只 在 第 6 个 脉冲 时 序 将 公共 
噪声 源 注入 到 系统 中 ,其 他 时 序 情况 下 不 注入 。 同 
理 , 内 部 噪声 源 配置 的 控制 开关 时 序 为 L000000 
1100]. Æ% 7,8 两 个 时 序 将 内 部 噪声 源 注 入 系 
统 , 其 他 情况 下 不 注入 。 仿 真 中 不 同 噪声 源 对 应 的 
迪克 开关 脉冲 时 序 设 置 情况 如 表 3 所 示 。 


表 3 系统 开关 时 序 控制 表 


Table 3 Sequential control of the system switch 


脉冲 开关 控制 时 序 
匹配 负载 0 0 0 0 0 0 0 1 1 1 
外 部 噪声 源 0 0 0 0 0 1 0 0 0 0 
内 部 噪声 源 0 0 0 0 0 0 1 1 0 0 
T Tz T 
Ta T5 
系统 噪声 源 HT FT, HT 
F T. F T ree 
TT "ee +T 


总 之 ,通过 上 述 多 个 开关 脉冲 时 序 的 搭配 设置 ， 
对 系统 注入 不 同 的 噪声 源 ,以 实现 期 望 的 几 种 输入 
噪声 状态 。 
5.1.2 两 点 定 标 处 理 

在 中 频 单元 的 耦合 输出 端 观测 全 功率 信和 号 的 波 
形 ,并 且 将 输出 值 换算 为 对 应 的 噪声 源 亮 温 值 。 在 
换算 过 程 中 ,涉及 到 辐射 计 系统 的 两 点 定 标 问题 。 
对 仿真 系统 注入 已 知 温度 的 噪声 源 及 匹配 负载 , 进 
行 两 点 定 标 ,获得 系统 的 亮 温 一 功率 定 标 方程 。 在 
仿真 过 程 中 ,系统 首先 注入 匹配 负载 噪声 源 ,再 考虑 
到 接收 机 本 身 的 噪声 ,获得 此 时 整个 系统 的 亮 温 值 
T, —500 天 ,对 测量 得 到 的 中 频 耦 合 输出 功率 值 求 
平均 ,得 到 功率 Pi 一 1. 6489 10 ”W; 接 着 在 以 上 
系统 中 添加 内 部 噪声 源 ,整个 系统 对 应 的 亮 温 值 改 
变 为 T: =600 天 ,中 频 耦 合 输出 功率 求 平 均 后 得 到 
P,—1.9803X10 “”W。 通过 以 上 两 组 数据 进行 两 
点 定 标 ,从 而 确定 接收 机 系统 的 亮 温 一 功率 定 标 线 ， 


接收 机 系统 的 定 标 线 如 图 10 所 示 。 
图 10 中 ,直线 的 斜率 对 应 系统 的 增益 值 ,表征 
接收 机 系统 亮 温 值 与 信号 输出 功率 之 间 的 线性 
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图 10 接收 机 系统 定 标 线 


Fig. 10 Calibration line of the receiver system 
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5.1.3 系统 指标 参数 分 析 

获得 定 标 线 以 后 ,对 于 未 知 噪声 亮 温 值 的 天 线 
信号 ,就 可 以 通过 测量 其 在 接收 机 系统 输出 端的 输 
出 功率 ,在 该 直线 上 进行 亮 温 值 的 换算 。 进 行 功率 
换算 以 后 的 中 频 耦 合 输出 的 单 通道 全 功率 信和 号 波形 
图 如 图 11 rz. 

在 图 11 中 ,上 半 部 分 表示 中 频 耦 合 输出 的 功率 
值 , 下 半 部 分 表示 换算 以 后 的 系统 亮 温 值 , 图 中 的 黑 
线 表 示 仿 真 输出 , 灰 线 表示 理论 输出 。 系 统 切 换 于 
5 种 噪声 源 注 人 状态 :第 一 种 注 和 天线 和 接收 机 只 
声 信号 ,等 效 总 亮 温 值 为 280 K; 第 二 种 注 和 人 天线 、 
公共 噪声 源 以 及 接收 机 噪声 信号 ,等 效 总 亮 温 值 为 
350 K; 第 三 种 注 和 人 天线、 内 部 噪声 源 以 及 接收 机 噪 
声 信和 号 ,等 效 总 亮 温 值 为 380 K; 第 四 种 注入 匹配 负 
载 、 内 部 噪声 源 以 及 接收 机 噪声 信号 ,等 效 总 亮 温 值 
为 600 开 ; 第 五 种 注入 匹配 负载 以 及 接收 机 噪声 信 
号 ,等 效 总 亮 温 值 为 500 K。 从 图 中 可 以 看 出 ,通过 
对 时 序 脉冲 的 控制 ,实现 了 不 同 噪声 源 比 例 的 注入 ， 


思 等 ;基于 Simulink 的 干涉 式微 波 辐射 计 系统 仿 直 研究 763 


将 系统 切换 于 不 同 的 状态 之 间 。 对 仿真 中 每 个 系统 
状态 下 输出 的 亮 温 值 ,经 过 求 平均 及 均 方差 值 处 理 ， 
分 别 与 理想 数据 进行 对 比 ,结果 如 表 4 所 示 。 


N 3 
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R 
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y 000 
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一 一 理论 输出 ”一 一 仿真 输出 


图 11 中 频 耦 合 输出 图 
Fig. 11 Coupling output of IF 


表 4 全 功率 测量 仿真 数据 对 比 


Table 4 Comparison of the simulation data 


数据 状态 1 状态 2 状态 3 状态 4 状态 5 
Tsys/K 280 350 380 600 500 
仿真 亮 温 均值 /K 278. 8597 348. 7614 378. 8509 597. 5534 497. 5705 
AT/K 8. 8544 11. 0680 12. 0167 18. 9737 15. 8114 
仿真 亮 温 均 方 差 /K 9. 1030 11. 2928 12. 1044 19. 1612 15. 8415 


表 4 中 ,理想 灵敏 度 AT 的 值 是 通过 式 (1) 进 行 
计算 所 得 ,在 仿真 系统 中 ,B 为 2 kHz,r 为 0.5s。 
通过 数据 对 比 可 以 看 出 仿真 输出 亮 温 平均 值 与 理论 
值 基本 一 致 ,相对 误差 在 0.5% 左 右 ; 均 方差 值 与 理 
论 灵敏 度 误差 在 2.5% 之 内 ,理论 值 与 仿真 值 相差 
不 大 ,验证 了 仿真 系统 的 正确 性 。 

需要 说 明 的 是 ,此 处 在 计算 灵敏 度 参数 时 使 用 
的 + 为 系统 的 单 点 积分 时 间 , 实 际 仿真 是 在 一 个 较 
长 的 观测 周期 内 对 单 点 积分 时 间 进行 统计 平均 , 通 
过 拉 长 系统 的 积分 时 间 , 达 到 优化 系统 灵敏 度 参数 
的 目的 。 仿 真 中 设 定 观测 周期 为 100 s, 将 单 点 积 4 
时 间 所 对 应 的 灵敏 度数 值 优 化 了 V200 倍 。 在 表 4 
对 比 过 程 中 ,为 了 重点 讨论 说 明 仿真 所 得 与 理论 灵 
敏 度 参数 的 对 应 关系 ,所 以 只 比较 了 单 点 积分 时 间 
所 对 应 的 灵敏 度数 值 。 

接 下 来 考察 系统 的 线性 度 指标 ,通过 仿真 系统 
的 输入 噪声 温度 值 与 输出 功率 值 之 间 的 线性 关系 进 
行 说 明 。 仿 真 系统 中 设置 了 几 组 不 同 的 亮 温 输入 


值 ,分别 得 到 了 对 应 的 功率 输出 值 , 利 用 所 得 数据 绘 
制 出 辐射 计 线性 度 仿真 曲线 如 图 12 所 示 。 
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Fig.12 Linearity of the silulation system 
由 图 12 可 以 看 出 ,仿真 模型 的 输出 值 与 输入 的 
噪声 温度 之 间 具 有 和 良好 的 线性 关系 ,其 线性 相关 度 
用 线性 相关 系数 尽 来 表示 : 


Ep V= LY V, 表示 输出 功率 的 平均 值 ,T = 


7 一 1 
LIT ,表示 对 应 噪声 温度 的 平均 值 
1 一 0 


将 图 12 中 的 数据 分 别 代 入 式 (2) 中 ,计算 得 到 
搭建 的 辐射 计 仿 真 系统 的 线性 相关 系数 为 R = 
0.9999, 结果 表明 仿真 系统 的 输出 值 与 输入 值 之 间 
具有 良好 的 线性 度 。 

5.2 双 通 道 干涉 测量 仿真 

对 于 双 通 道 干 涉 测量 ,首先 简化 系统 注入 噪声 
源 的 状态 ,只 把 天 线 信号 作为 公共 信号 注 人 到 两 个 
接收 机 通道 中 , 记 为 Ti ,同时 ,接收 机 自身 的 噪声 
言 号 作为 各 个 通道 的 系统 噪声 源 ,分 别 记 为 Tua 和 
Te。 在 理论 上 ,系统 输出 的 理论 模拟 相关 系数 的 幅 
值 大 小 与 噪声 源 的 分 配 值 有 如 下 关系 : 

T uie 
dee Lea o daa T Tag 

通过 在 系统 中 配置 不 同比 例 的 噪声 源 , 实 现 模 
拟 相 关系 数 不 同 幅 值 的 输出 . 另 一 方面 ,对 于 双 通道 
的 二 元 干涉 仪 仿真 系统 ,还 涉及 到 两 个 接收 机 通道 
之 间 的 相位 差 问题 ,这 个 差 值 与 得 到 的 相关 系数 相 
位 值 是 一 一 对 应 关系 。 

利用 仿真 系统 输出 数字 相关 系数 7, 经 过 对 干 
涉 测量 输出 值 的 校正 ,实现 数字 相关 系数 到 模拟 相 
关系 数 的 估计 ,得 到 模拟 相关 系数 osim ,对比 不 同 


(3) 


Pideal 


中 ,主要 从 均值 和 均 方差 两 个 角度 进行 结果 评估 ,一 
方面 ,仿真 输出 的 模拟 相关 系数 在 幅度 和 相位 上 的 
均值 应 该 与 理论 值 相近 ; 男 一 方面 ,理论 上 亮 温 值 的 
均 方差 对 应 着 系统 灵敏 度 的 大 小 , 则 利用 式 (3) 可 
以 得 到 模拟 相关 系数 均 方 差 值 的 理论 值 应 该 是 
1/ VBr ,对 比 此 理论 值 与 实际 仿真 输出 的 模拟 相 
关系 数 的 均 方差 值 ,评估 仿真 性 能 。 
5.2.1 相关 系数 仿真 评估 

将 两 个 接收 机 通道 的 相位 差 值 固定 为 0 ,接收 机 
通道 噪声 源 Ta 和 Tree 固定 为 300 KK, 只 改变 天 线 噪声 
信号 Tus 值 大 小 ,通过 仿真 输出 测量 oaa 与 oss 数据 。 


示例 ,设置 天 线 噪声 信号 亮 温 为 Tue = 75 K, 
根据 式 (3), 有 [oua | = 0.2。 系 统 输出 数字 相关 系 


数 ,再 利用 数字 相关 系数 到 模拟 相关 系数 的 估计 理 
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论 ,得 到 估计 的 模拟 相关 系数 值 , 估 计 得 到 的 模拟 相 
关系 数 如 图 13 所 示 。 


0.4 


0 20 40 60 80 100 120 140 160 180 200 
时 间 /s 
(a) 估计 得 到 的 模拟 相关 系数 幅 值 


20 40 60 80 100 120 140 160 180 200 
时 间 /s 
(b) 估计 得 到 的 模拟 相关 系数 相位 


图 13 估计 得 到 的 模拟 相关 系数 图 


Fig. 13 Estimated analog correlation coefficient 


对 所 得 的 仿真 数据 结果 进行 处 理 分 析 ,分 别 求 
得 模拟 相关 系数 幅 值 和 相位 的 均值 与 均 方 根 值 ,得 
到 仿真 输出 的 模拟 相关 系数 | oss. | 均值 为 0. 2033. 
而 理想 模拟 相关 系数 | oua | 为 0.2, 相 对 误差 为 
1. 65 76 ;模拟 相关 系数 /om 均值 为 一 0. 2932 ,理想 
模拟 相关 系数 一 ou 为 0 。 差 值 控制 在 0. 3 以 内 ; 
仿真 相关 系数 oum 均 方 根 值 为 0.0380 ,理论 模拟 相 
关系 数 ou 的 均 方 根 值 为 0.0322, 两 者 较为 接近 。 
通过 数值 对 比 看 出 仿真 输出 与 理论 值 差 异 不 大 ,说 
明了 仿真 系统 的 正确 性 。 

实验 仿真 了 多 组 不 同 分 配 比 例 的 Ta 和 Tree， 
分 别 求 得 对 应 的 均值 及 均 方差 值 ,并 与 理论 值 进 行 
对 比 , 结 果 如 表 5 Bro. 

RS 改变 噪声 亮 温 值 结果 对 比 
Table 5 Result of different bright temperature 


数据 设置 1 设置 2 设置 3 设置 4 设置 5 
区 33. 3 75 100 200 300 
| piaea | 0.1 0.2 0. 25 0.4 0.5 
| osim | 0.0998 — 0.2033 0.2595 0.4431 0.5802 
Z pida / 0 
Z pin / ^ 0.1670 —0.2932 —0.3680 —0.2040 —0.1440 
iae BIJT 258 0. 0322 
psim 均 方差 0.0360 0.0380 0.0382 0.0392 0.0398 


通过 表 5 可 以 看 出 ,仿真 得 到 的 模拟 相关 系数 
pin 与 理论 模拟 相关 系数 cum 幅度 均值 相对 误差 在 
5% 左 右 , 相 位 误差 控制 在 0. 企 以 内 , 均 方 差 值 也 具有 
很 好 的 一 致 性 。 对 比 结果 表明 仿真 过 程 的 正确 性 。 
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5.2.2 相关 相位 仿真 评估 
国定 Tus Te 大 小 比例 的 分 配 , 只 对 两 个 接收 
机 通道 的 相位 差 值 进行 改变 ,观测 仿真 输出 的 oua 
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与 Psim 的 差异 。 Tua 7775 K, Tra —Tu — 300 K, 则 
理论 模拟 相关 系数 Lou] 29 0.2. TFIEESETI T fH 
仿真 实验 ,得 到 的 数据 结果 如 表 6 所 示 。 


表 6 改变 接收 机 相位 差 结果 对 比 


Table 6 Result of different phase of reciever 


数据 设置 1 设置 2 设置 3 设置 4 设置 5 设置 6 设置 7 
| Oideal | 0.2 
| psin | 0. 2033 0. 2029 0. 2034 0. 2033 0. 2029 0. 2034 0. 2033 
Z pica /^ 0 30 60 90 120 150 180 
Z sim /" — 0. 2932 29. 6953 59. 7082 89. 7068 119.6953 149.70182 179. 7068 
pideal 均 方差 0.0322 
psim 均 方差 0. 0380 0. 0372 0. 0377 0. 0380 0. 0372 0. 0377 0. 0380 


通过 表 6 数据 可 以 看 出 ,仿真 得 到 的 模拟 相关 
系数 ps 与 理论 模拟 相关 系数 oum 幅 度 均值 相对 误 
差 为 1.7% 以 内 ,相位 误差 控制 在 0. 3 左右 , 均 方 差 
值 也 很 接近 。 综 合 对 比 结果 表明 了 仿真 的 正确 性 。 
5.3 ADC 采样 频率 对 相关 结果 的 影响 

对 干涉 式微 波 辐射 计 系统 来 说 ,ADC 部 分 对 信 
号 的 处 理 是 得 到 复 相关 输出 的 关键 环节 。 在 ADC 
中 ,涉及 到 对 原始 信号 的 采样 问题 ,要 想 获得 信号 完 
整 的 重建 ,采样 时 钟 频率 必须 满足 奈 奎 斯 特 准则 。 
对 于 所 搭建 的 仿真 系统 ,采用 了 正 交 混 频 方法 ,将 中 
频 输出 信和 号 分 为 /Q 两 路 , 则 理论 上 只 需要 一 半 的 
奈奈 斯 特 采 样 率 即 可 保证 不 发 生 信 号 的 频 域 混 炙 。 
理论 上 ,不 同 的 时 钟 采样 频率 对 相关 结果 会 产生 影 
Ul] ,采样 频率 越 高 ,信号 的 重建 越 准 确 , 输 出 的 相关 
结果 指标 会 越 好 。 

仿真 系统 的 中 频带 宽 为 2 kHz, 则 可 以 确定 对 
信号 进行 完整 重 现 的 最 低 采样 频率 是 2 kKHz。 固 定 
天 线 噪声 信号 的 亮 温 值 为 Ti 二 75 K ,接收 机 噪声 
记号 的 亮 温 值 为 Trea — T««—300 开 , 接 收 机 之 间 相 
位 差 为 0 ,对 仿真 系统 设置 不 同 的 采样 频率 人 , 测 
量 相关 输出 的 均值 及 均 方 差 值 。 

数据 结果 显示 , 随 着 采样 率 的 变化 ,相关 系数 的 
均值 变化 不 大 ,但 均 方 差 值 有 所 影响 。 利 用 所 得 数 
据 绘制 图 形 ,观测 采样 率 对 相关 系数 均 方 差 值 的 影 
响 趋势 , 横 轴 取 采 样 率 与 系统 带宽 的 比值 , 即 人 /B， 
纵 轴 取 理 论 均 方差 与 仿真 输出 均 方差 的 比值 , 即 
Qiaca /psim， 观 察 仿真 输出 数据 随 采 样 率 的 变化 ,得 到 


的 图 形 如 图 


图 14 中 显示 的 散 点 就 是 仿真 获得 的 数据 , 通 


14 所 示 。 


据 拟 合 以 后 产生 图 中 所 示 曲 线 。 由 上 图 可 以 看 出 , 随 
着 系统 采样 率 的 增加 ,仿真 输出 相关 值 的 均 方差 


越 趋 近 于 理论 值 


越 来 


,表明 系统 的 灵敏 度 指标 不 断 提升 。 


i 


同时 ,曲线 前 期 上 升 趋势 比较 明显 ,而 后 期 趋 近 于 3 
缓 ,说 明 在 满足 奈 夺 斯 特 准 则 的 条 件 下 ,初始 采样 率 的 
增加 对 系统 性 能 有 明显 提升 效果 ,而 当 采样 率 达到 一 
定 指标 之 后 ,改善 系统 灵敏 度 的 效果 不 再 显著 。 


"0 1 2 3 4 5 6 7 8 
采样 率 /带宽 


* 仿真 数据 数据 拟 合 
E] 14 采样 率 对 相关 输出 的 影响 图 
Fig. 14 The influence of sampling rate on the correlation 
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本 文 所 做 的 主要 工作 是 利用 Matlab/Simulink 
平台 搭建 出 干涉 式微 波 辐 射 计 系统 ,对 二 元 干涉 仪 
进行 了 时 域 上 的 仿真 。 利 用 所 搭建 的 模型 实现 了 单 
通道 全 功率 测量 仿真 以 及 双 通 道 干涉 测量 仿真 , 完 
成 了 系统 定 标 三 阶 量化 数字 相关 结果 的 转化 以 及 
ADC 采样 频率 影响 的 评估 ,在 时 域 范围 内 实现 了 对 
系统 测量 精度 .灵敏度 指标 的 定量 化 分 析 。 
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Simulation of Microwave Interferometric Radiometer 


System based on Simulink 
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Abstract: A method using Simulink to build the simulation model of the microwave interferometric radiome- 


ter (MIR) system is introduced in this paper. T'wo-element interferometer is the basic unit of the MIR sys- 


tem. By the Simulink platform. a simulation model of the two-element interferometer was built and the 


measurement process of the random noise was simulated. Through the simulation of the full power meas- 


urement achieved interference measurement of single channel and two channels. Quantitative analysis of the 


radiometric accuracy and radiometric sensitivity were achieved in the time domain. 


Key words: Microwave interferometric radiometer; Visibility function; Simulation model; Simulink 


